High-frequency flow reversal of AC electro-osmosis due to steric effects 
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The current theory of alternating-current electro-osmosis (ACEO) is unable to explain the experi- 
mentally observed flow reversal of planar ACEO pumps at high frequency (above the peak, typically 
fO — fOO kHz), low salt concentration (1 — fOOO (jM), and moderate voltage (2 — 6 V), even if taking 
into account Faradaic surface reactions, nonlinear double-layer capacitance and bulk electrothermal 
flows. We attribute this failure to the breakdown of the classical Poisson-Boltzmann model of the 
diffuse double layer, which assumes a dilute solution of point-like ions. In spite of low bulk salt 
concentration, the large voltage induced across the double layer leads to crowding of the ions and 
a related decrease in surface capacitance. Using several mean-field models for finite-sized ions, we 
show that steric effects generally lead to high frequency flow reversal of ACEO pumps, similar to 
experiments. For quantitative agreement, however, an unrealistically large effective ion size (several 
nm) must be used, which we attribute to neglected correlation effects. 



I. INTRODUCTION 
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Electrokinetic phenomena, which couple fluid flow, ion 
transport and electric fields in electrolytes, are exploited 
in a variety of microfiuidic technologies A solid sur- 
face in contact with an electrolyte typically acquires a 
surface charge and forms an electric double layer com- 
posed of wall charge and a screening layer of excess coun- 
tcrions. Electro-osmosis results from the action of an ap- 
plied electric field on this double-layer charge, which sets 
the bulk fluid into motion. Although the double layer is 
typically very thin (nm) compared to the channel dimen- 
sions (jum), it is the source of long-range flow patterns. 

The classical continuum theory of electro-osmosis 0, H[ 
suffices to understand most direct-current (DC) electro- 
osmotic flows, which are linear in the applied voltage. 
DC electroosmosis results from application of an external 
electric field parallel to the surface, which acts on equi- 
librium double-layer charge and induces motion in the 
bulk, neutral liquid. The linear response is due to the 
assumption that the applied field does not perturb the 
pre-existing surface charge. When DC electroosmosis is 
used in microfiuidic environments, electrodes are placed 
at the ends of relatively long (cm) channels, and there- 
fore large voltages (kV) are typically needed to generate 
sufficient fields. 

Electrokinetic phenomena can be exploited at much 
lower voltages and with much greater precision by plac- 
ing electrodes close together inside the fluid channels and 
using alternating current (AC) to inhibit Faradaic reac- 
tions. Large electric fields (100 V/cm) needed to induce 
adequate flow can thus be generated with only a few 
volts by micro-electrodes. The most studied example is 
alternating-current electro-osmosis (ACEO) [4|, discov- 
ered by Ramos et al. [E IE 0, @] ■ When a periodic array 
of interdigitated electrodes is placed inside the channel, 
an AC signal applied to the electrodes can generate a 
steady set of micro-vortices. Ajdari showed that if ge- 
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FIG. 1: One period of an asymmetric planar array of 
microelectrodes in an AC electro-osmotic (ACEO) pump. 
This geometry has been studied experimentally by several 
groups [iol [la . [l4l . [l8T ] and is the subject of the present 
theoretical study. The dimensions are; Wl = 4.2 /im , 
W2 = 25.7 /im, Gl = 4.5 urn , and G*2 = 15.6 /im, using 
the notation of Ref. [l9| . 



ometrical asymmetry is introduced, then these vortices 
can be rectified to drive net pumping in one direction 
over the array [9j. The original implementation of this 
principle was a planar array of flat electrode pairs of un- 
equal widths and gaps [13, [ill, EH' as shown in Figure 
[TJ Recent designs for "3D ACEO" pumps have achieved 
much faster flows (mm/sec) using non-planar electrodes, 
which more efficiently rectify opposing slip velocities in 
a "fluid conveyor belt" [HI, 0, [lH . Waves of voltage can 
also pump fluids over electrode arrays by traveling- wave 
electro-osmosis (TWEO) [MlIS- 

These effects exemplify the fundamental nonlinear 
electrokinetic phenomeno n 12011 of induced-charge electro- 
osmosis (ICEO) [U, [22|, |23|. which also has other ap- 
plications in microfluidics. The key source of nonlin- 
earity is that the diffuse double-layer charge is induced 
near a polarizable surface by the applied electric field, 
which then acts on it to drive nonlinear electro-osmotic 
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flow. Although ICEO flows were first described in the 
1980s in the context of colloid science 0, [Hj], they are 
finding new applications in microfluidic devices. As in 
the case of ACEO pumps, broken spatial symmetries in 
ICEO flows can also be used to pump fluids around po- 
larizable microstructures and to manipulate polarizable 

particles [Hm. 

The current theory of ICEO flow [13], including 

ACEO @, [H and TWEO [H, E3 microfluidic pumps, 
is based on the classical electrokinetic equations @, yj]- 
This century-old model comprises the Poisson-Nernst- 
Planck (PNP) equations of ion transport coupled to the 
Navier-Stokes equations of viscous fluid flow via elec- 
trostatic Maxwell stresses. The crucial modification is 
in the boundary conditions, relaxing the assumption of 
fixed surface charge (or constant zeta potential) to allow 
for significant electrostatic polarizability of the surface 
(e.g. fixed surface potential). The electrokinetic equa- 
tions themselves, however, have not been questioned, un- 
til very recently (29| . which provides the motivation for 
this work. 

ICEO flows typically involve large voltages induced 
across the double layer, greatly exceeding the thermal 
equilibrium voltage, kT/e (= 25 mV at room tempera- 
ture). In particular, experimental ACEO electro-array 
pumps [IS El HEED!], the subject of our study, usu- 
ally apply several volts to the double layer. Such large 
voltages of order lOOfcT/e inevitably lead to the break- 
down of the PNP equations by violating the fundamental 
assumption of an ideal dilute solution [30L l3lj . The limi- 
tations of dilute solution theory for double-layer structure 
are well known, and many attempts have been made to 
incorporate effects such as steric exclusion and electro- 
static correlations. (See Ref. [3(| for a recent review.) 
ICEO flow, however, raises new issues due to its more 
extreme, dynamical context 29]. Even if the bulk solu- 
tion is dilute, the surface becomes so highly charged by 
the induced voltage that counterions become crowded in 
the double layer, while simultaneously driving tangential 
fluid flow. 

The need for improved models is also indicated by some 
unexplained features of ICEO experiments. In this paper, 
we focus on the tendency for planar ACEO pumps to re- 
verse at high voltage [H, El, El, HI, [HI , which naturally 
poses problems for microfluidic applications, if it can- 
not be reliably predicted. The standard model of ACEO 
based on the classical electrokinetic equations predicts a 
single peak in the flow rate versus frequency at the char- 
acteristic "RC" charging frequency of the electrodes, in 
both the linear Q and nonlinear [H (34[ regimes. Flow 
reversal has been observed at high voltage (> 2 V) and 
high frequency (10-100 kHz) in ACEO pumping of dilute 
KC1 fll^.fl8l.l32| and deionized water [3] with 10/zm scale 
electrode arrays. 

Representative examples of the observed and predicted 
frequency response are shown in Figure [2] for the pump 
geometry of Figure [1] Here we compare experimental 
data taken at 3.0 V = 120fcT/e in an 0.03 mM KC1 so- 
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FIG. 2: Pumping velocity vs. frequency for the planar ACEO 
pump in Figure [1] from experiments and simulations. Points: 
Experimental data at 3 Vpp in 0.03 mM KC1 from Ref. [Hj] 
(taken by J. P. Urbanski); Solid curve: Linear Debye-Hiickel- 
Stern model of the diffuse layer with a constant compact layer 
capacitance C s = e/(5\d), where 5 — 0.5 is chosen to fit 
the experimental peak frequency. Dashed curve: Nonlin- 
ear Gouy-Chapman-Stern (Poisson-Boltzmann) model of the 
diffuse layer also with 5 — 0.5. 

lution over gold electrodes [jl| to simulations with two 
commonly used models: one assumes the linear Debye- 
Huckel (DH) model of the double layer 0, 0, H, while 
the other adopts a non-linear Poisson-Boltzmann (PB) 
model based on the Gouy-Chapman solution [19]. Fol- 
lowing prior work [1,(1 [231] , a constant Stern-layer capaci- 
tance is added in order to better fit the experimental data 
by matching the peak frequency in the DH model. Both 
models can reasonably fit the data in terms of the peak 
in the forward pumping, though this becomes more diffi- 
cult for PB at higher voltages [l9j . Each model predicts 
a single peak of forward flow, while the experiments show 
the pump reversing above 10 kHz. The inability of these 
classical models to capture, even qualitatively, this high 
frequency flow reversal motivates us to consider modified 
double- layer models. 

The reversal of ACEO pumps has been considered 
in several previous studies. Flow reversal was first at- 
tributed to Faradaic reactions [32], and this view of 
"Faradaic charging" persists in recent work on DC- 
biased AC electro-osmosis (35|. Simulations of ACEO 
pumps with Butler- Volmer reaction kinetics, however, 
have failed to predict the observed flow, especially at 
high frequency [19l. [33], In existing models, a weak flow 
reversal due to reactions can only be observed at low 
frequency (far below the maximum) and for certain sets 
of parameters [1 El [HI- This is consistent with the 
observation of gas bubbles from electrolysis at low fre- 
quency and high voltage in dilute KC1 with gold elec- 
trodes jTij. Recently, higher resolution measurements of 
the same pump design with platinum electrodes have re- 
vealed very weak (< 10/im/sec) reverse ACEO flow at 
low frequency (< 20 kHz) and low voltage (< 1.5 V) and 
have demonstrated the importance of Faradaic currents 
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through in situ impedance spectroscopy [36|. Although 
the theory predicts similar flow reversal, it is not in quan- 
titative agreement and does not predict the concentration 
dependence. 

In addition to Faradaic reactions, various other non- 
linear effects in dilute solution theory also dominate at 
low frequencies: The differential capacitance of the dif- 
fuse layer diverges, which causes the RC charging time 
to grow exponentially with voltage [HI , and salt adsorp- 
tion and tangential conduction by the diffuse layer are 
coupled to (much slower) bulk diffusion [3?], [M, [39[ . All 
of these effects have recently been incorporated in sim- 
ulations of ACEO at large voltages using the classical 
electrokinetic equations, but high-frequency flow rever- 
sal was not observed ■ Another pos sible source 
of flow reversal is AC electrothermal flow |4Cj , which can 
lead to reverse pumping in experiments with planar elec- 
trode arrays [35|], but under different conditions of much 
higher salt concentration (> 1 M), voltage (> 10 V) and 
higher frequency (> 100 kHz) than ACEO (< 0.01 M, 
< 5 V, < 100 kHz). To date, none of the above ef- 
fects have been able to explain the experimental data for 
ACEO pumps. 

In this paper, we propose steric effects of ion crowding 
in the double layer as a possible cause of the observed 
flow reversal at high frequency. Using simple modifica- 
tions of the electrokinetic equations to account for finite 
ion sizes, we are able to predict a frequency response 
of ACEO pumps that is very similar to that observed 
in experiments. Although our results do not provide a 
complete quantitative theory, we will demonstrate that 
accounting for steric effects can have an important im- 
pact on the theory of ACEO. 



II. STERIC EFFECTS IN A THIN DOUBLE 
LAYER 

In the dilute-solution theory of electrolytes @, EJ , the 
concentration of each ionic species in the diffuse part of 
the electric double layer at a charged surface is in thermal 
equilibrium with a Boltzmann distribution, 



Ci(x) = c exp 



-z t e4>(x) 
kT 



(1) 



where k is Boltzmann's constant, T the temperature, e 
the electron charge, Zi the valence, cq the concentration in 
the neutral bulk electrolyte just outside the double layer, 
and 4> is the (mean) electrostatic potential relative to the 
bulk (x S> \d, where Ad is the Debye-Hiickel screening 
length) . In classical linear electrokinetic phenomena, the 
diffuse-layer voltage drop 4>(0) is set by chemical equi- 
librium at the surface and is thus typically comparable 
to the thermal voltage kT/e. At an electrode driving 
ACEO flow, however, much larger diffuse-layer voltages 
(f)(0) S> kT/e are induced by a applied voltage of order 
lOOfcT/e. Under these conditions, it is easy to see that 



Boltzmann equilibrium ([T]) breaks down by predicting di- 
verging concentrations of countcrions at the surface, of 
order e 100 co. 

The assumption of a dilute solution is thus incompat- 
ible with a large applied voltage. If we assume ions have 
a characteristic length scale, a, then the corresponding 
cutoff concentration a~ 3 for the breakdown of dilute- 
solution theory is reached at relatively low voltage, even 
if the bulk salt solution is very dilute. For example, for 
Co = 10 -5 M, z = 1, and a = 3 A(including a solvation 
shell), the cutoff concentration is reached at 0.33 V. To 
account for the excess ions (at typical equilibrium volt- 
ages), Stern postulated a compact layer of solvated ions 
of finite size on the surface [4l| . which carries most of 
the double-layer voltage as the diffuse-layer capacitance 
diverges. Such an intrinsic "surface capacitance" is also 
invoked in models of ICEO flows, where it can also in- 
clude the effect of a thin dielectric coating on a metal 
surface 0, It seems unlikely, however, that an atom- 
ically thin Stern or coating layer could withstand sev- 
eral Volts, e.g. since dielectric breakdown occurs in most 
materials (including water) in fields of order 10 MV/m 
= 0.01 V/nm [42j ■ Instead, it seems the diffuse layer must 
carry a substantial voltage <j> ^> kT/e in ACEO experi- 
ments, which causes the ions to become highly crowded 
near the surface. 

A variety of "modified Poisson-Boltzmann equations" 
(MPB) have been proposed to describe equilibrium ion 
profiles near a charged wall, as reviewed in Ref. [30j . To 
capture ion crowding effects across a wide range of volt- 
ages, we employ the simplest possible MPB model [13, EH, 
EL EI| ) first proposed by Bikerman [46| , which is a con- 
tinuum approximation of the entropy of ions on a lattice 
of size a. As shown in Fig. a), when a large voltage is 
applied, the counterion concentration exhibits a smooth 
transition from an outer PB profile to a condensed layer 
at c = c max = a~ 3 near the surface. 

In ACEO, charging dynamics are very important, so 
the double-layer capacitance is an important property. 
Dilute-solution theory predicts that the differential ca- 
pacitance, Cd, diverges with the voltage as, 
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where vPd is the voltage applied across the double layer. 
For a concentrated solution described by Bikerman's 
model, we predict the opposite trend [30(, 



X77 sinh ( 



kT I 



1v sinh 



, 2kT , 



|ln[l + 2^sinh 2 (^)] 
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where v = 2coa 3 is the bulk volume fraction of ions [3C 
(Note that the same formula was derived by Korny- 
shev [47j in the context of ionic liquids around the same 
time that Kilic et al. [3(1 derived it in the present context 
of electrolytes.) 

As shown in Fig. [3jb), with this model the capaci- 
tance reaches a maximum near the critical voltage then 
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FIG. 3: (a) Schematic of the equilibrium distribution of coun- 
terions near a positively charged surface taking into account 
a minimum ion spacing a for large applied voltages. As the 
voltage increases, the width of the condensed layer will in- 
crease causing a decrease the capacitance, (b) The voltage 
dependence of the capacitance of the diffuse layer for PB and 
MPB at different values of v\ where v = 2a 3 co is the bulk 
volume fraction of ions. 



decreases at larger voltages because the capacitor width 
grows due to steric effects. The same result also fol- 
lows from a simpler model using PB theory outside a 
condensed layer of variable width at the maximum den- 
sity [30|| . We will show that this general decrease in 
capacitance with voltage that can drastically impact the 
charging dynamics in ACEO flows. 

As Biesheuvel et al. point out [HI, H^], Bikerman's 
latticed-based approach can underestimate the true steric 
effect to a large degree, and more accurate models of the 
entropy of a hard-sphere liquid, such as the Carnahan- 
Starling (CS) equation of state or various extensions to 
multicomponent mixtures [49|, [5(|, can be used instead. 
The CS equation is not tractable for a closed-form analyt- 
ical expression for double-layer capacitance, but we can 
easily obtain it numerically as a function of the voltage. 
For the details, see the Appendix. 

Although various mean-field models, such as CS, are 
very accurate for equilibrium properties of the hard- 
sphere gas and liquid at low to moderate volume frac- 
tions, $ < 0.55, none can claim to accurately describe 
the packing limit at large volume fractions, which would 



generally lead to the decrease in double-layer capacitance 
at large voltage [30| . For example, these models predict 
a diverging chemical potential only at $ = 1, and not 
and not at the mathematical packing bound (<& = 0.74, 
corresponding to the face-centered cubic lattice) or the 
more relevant bound of the jamming transition (around 
$ = 0.63 for hard spheres [5l|). Even in simple model 
systems in equilibrium, the jamming transition is poorly 
understood [52|, HU, Q, l55l . [5|| , so none of these mod- 
els can possibly be expected to accurately capture the 
dynamics of jamming of hydrated ions very close to a 
surface driven by a large AC voltage, while also produc- 
ing ACEO flow. Of course, under such extreme condi- 
tions, the mean-field approximation for both steric and 
electrostatic interactions should also be questioned. 

Due to all of these complexities, following Bazant et 
al. [2^ |. we will focus on the simplest approximation of 
Bikerman's model in order to explore the qualitative im- 
pact of steric effects on ACEO flow. It will turn out 
that this very simple model suffices to predict high fre- 
quency flow reversal of ACEO pumps. We will then re- 
peat some calculations with the CS model, which im- 
proves the agreement with experiments somewhat, but 
does not alter the qualitative predictions of Bikerman's 
model. 



III. MODEL FOR ACEO PUMPS 

In our theoretical study of ACEO pumping, we con- 
sider the standard experimental geometry of planar, 
asymmetric electrode-array pumps shown in Fig. [T] due 
to Brown et al. 10]. The specific dimensions are selected 
to match the experiments Studer et. al [l2| and Urban- 
ski et al. [3, Il8j| . The basic question of ACEO flow 
in such pumps has been thoroughly addressed by Olesen 
et. al [191 l34j using the classical electrokinetic equations 
of dilute solution theory 0, [1[ . The significant difference 
between this work and Ref. [p| is our application of 
the new physical model for the double-layer capacitance, 
taking into account steric effects at large voltage. Unlike 
Ref. [3J] , we also ignore "strongly nonlinear" effects asso- 
ciated with salt uptake by a hi ghly charged diffuse layer 
and resulting bulk diffusion [37lJ3q] , since these phenom- 
ena are greatly reduced by steric effects [3(1 Hl| ■ 

In the bulk fluid outside the electric double layers, 
we assume electroneutrality with a constant electrical 
conductivity, a. Under these assumptions, the elec- 
tric potential in the fluid bulk, (j>, satisfies Laplace's 
equation, which is equivalent to Ohm's law for a con- 
stant bulk resistance. This corresponds to the "weakly 
nonlinear" regime where nonlinear circuit models still 
hold [m, [37], [38J . For boundary conditions on the poten- 
tial, we assume that there is no normal Ohmic current on 
the sections of channel substrate at y = and the upper 
boundary at y = H . 

For the boundary condition over the electrodes, we 
must incorporate a model of the thin electric double lay- 
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ers that will form. These double layers are so thin (10-100 
nm), that when we apply the boundary condition on <fi 
at y = 0, we assume this is outside the double layer and 
not at the electrode surface. The local charge stored in 
the double layer, q(x), changes in time due to the Ohmic 
current. Our double layer model provides the voltage 
difference between the fluid bulk and the externally ap- 
plied voltage at the electrode, V ex t, based on the amount 
of charge stored. Thus the boundary condition for the 
electric potential in the bulk is dynamically determined. 
While Faradaic current certainly play a role in ACEO, 
we neglect it in our model to focus on steric effects alone, 
which have not previously been considered. The role of 
reactions in classical double-layer models in ACEO was 
studied by Olesen et. al [l!| and shown not to play a role 
in high-frequency flow reversal. 

As in prior work 0, [l||, we introduce the following 
scales for the nondimensionalization of the governing 
equations, 

[x,y]=L, [t]=^, [fl = ^, [q] = —C . 

a ze ze 

Note that [t] is the "RC time" for the equivalent cir- 
cuit [H, [3?| , where the characteristic Debye capacitance 
is C = where e is the permittivity of the electrolyte 
and Xd is the Debye length. The characteristic length 
scale for computing the bulk resistance, L/a, is the elec- 
trode spacing L — G\. We assume no Stern layer with 
our new double layer model, though its inclusion would 
be trivial and would not significantly alter our results. 

The dimensionless problem reduces to Laplace's equa- 
tion for the potential in the fluid bulk, 



V 2 = 0. 



(4) 



The boundary condition on the substrate (the entire up- 
per wall at y — H/L and the insulating regions along 
y = between the electrodes) is 



dj)_ 
dy 



0, 



(5) 



The boundary condition over the electrode double layer 
is given by, 



(6) 



where ^D(q) is the total potential drop across the double 
layer. The functional relationship between the double 
layer voltage drop and charge is given by the model from 
Kilic et al. fH, 



(7) 



while the local charge is dynamically determined by the 
"RC" condition, 




dq d<fi 
dt dy 



The closed set of equations ((!])- © form our electrical 
model. Note that in the limit of very small v in Eq. 
[7J we obtain the classic result from Poisson-Boltzmann, 
^ d = — 2 sinh (I) . In the limit of small potential and 
small v we obtain ^ £> — —q, the linear Debye-Huckel 
model. The later two limits have been extensively studied 
[l9j : it is the impact of EqUJon the flow that is the focus 
of this paper. 

Once the electrical problem is solved, we compute the 
electroosmotic slip by the Helmholtz-Smoluchowski rela- 
tionship, 



U= 



94>{y = 0) 
dx 



(9) 



(8) 



and ignore any modifications at large voltage, e.g. due 
to crowding effects 29] , in order to focus on steric effects 
in the double-layer charging dynamics. Velocity is made 
dimensionless by the natural scale H [2l[ , 

eV 2 

where V is the magnitude of the applied voltage and r\ 
is the viscosity. An interesting fact, exploited by previ- 
ous authors @, H, 0], is that we can compute the net 
pumping without actually solving the flow field in the 
channel, but by taking the time and spatial average of 
equation [9] along the boundary y = 0. Further details of 
this derivation can be found elsewhere, e.g. in Olesen et. 
al [H|. 

Due to the irregular electrode geometry and the non- 
linear charging processes, we must resort to solving Equa- 
tions (j!])-© numerically. We couple the solution of 
Laplace's equation to the dynamic boundary condition 
for the charge. The solution is integrated forward in 
time until an adequate steady state is reached and the 
net flow is computed by averaging Eq. [5] in time and 
space. Details on the numerical method can be found in 
the supplementary material [59]. The simulations were 
confirmed by an independent method using the commer- 
cial software package COMSOL Multiphysics. 



IV. RESULTS 

We start by comparing our results with the steric 
model for the double layer, Eq. ([7]), to simulations of 
the dilute-solution models used in all prior work. In 
Fig. [4l we compare the flow velocity as a function of 
frequency for the classic linear Debye- Huckel model valid 
only when the applied voltage is much less that kT/e, 
the Poisson-Boltzmann non-linear capacitance which is 
recovered from our model when v — 0, and the model 
which accounts for steric effects, Eq. [7J At a voltage of 
lOOkT/e = 2.5 V, we see that the model which accounts 
for steric effect shows a high frequency flow reversal. As 
we note in Ref. L 29f] , to our knowledge, these results rep- 
resent the first demonstration of a mathematical model 
that can predict such reversal of flow in ACEO pumps. 
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FIG. 4: Average slip velocity as a function of frequency for 
three different models at 2.5 Volts. The solid blue curve is 
for v = 0.01 with the complete model, the black dash-dot line 
assumes the linear Debye-Huckel relationship q — —^d, the 
red dashed curve is for v = with a Stern layer capacitance 
of 8 = 0.1, Only the model with steric effects predicts flow in 
the negative direction. Points A, B, and C will be referred to 
in Fig. El 



In order to understand the physical mechanism of flow 
reversal, we must first understand why there is flow at 
all. The basic mechanism for flow in ACEO can be best 
understood when considering a geometrically symmet- 
ric electrode array, as in FigO subjected to a suddenly 
applied DC voltage (without any Faradaic current) [J]. 
When the voltage is initially switched on, the electric field 
lines point from the positive to the negative electrodes 
as shown in Fig[5]Ja). Since there is initially no dou- 
ble layer charge, the field lines are perpendicular at the 
electrode surface. A current will begin to flow through 
the electrolyte with ions migrating along the electric field 
lines and the double layers over the electrode will be- 
gin to charge as shown in Fig[5{b). The electric field is 
strongest at the electrode edges and the resulting high 
current quickly charges the double layer at that location. 
On the order of the characteristic charging time, t c , (as 
computed using the gap length scale) the electrode edges 
will be significantly screened. In Fig [3(b) we show the 
charge at y = at three instances in time. 

As the double layers charge up, the electric field lines 
in the fluid bulk are diverted around the screened elec- 
trode edges, resulting in a tangential component of the 
electric field. This is shown in Fig. [5jc) where we can 
see the field lines are no longer perpendicular over the 
entire electrode, especially near the electrode edges. The 
tangential electric field exerts a force on the fluid in the 
charged double layer, causing flow. Given the direction 
of the field and sign of the charge, the flow always goes 
from the edge of the electrodes, inward. The flow stream- 
lines and flow direction are shown in Fig. [5] (d). If the 
applied voltage is maintained, eventually the entire elec- 
trode becomes fully screened there is no electric field in 



the bulk and the flow will stop. The above argument 
about the direction of flow did not depend upon the sign 
of the applied voltage, reversing the sign causes the same 
flow pattern. Therefore, if an AC voltage is applied at a 
frequency corresponding to the charging time, time aver- 
aged vortex flow can be maintained. 

Ajdari introduced the idea of rectifying this symmet- 
ric flow pattern to create an ACEO pump by introducing 
various forms of geometrical asymmetry within each spa- 
tial period of the electrode array [9(. This principle was 
first implemented by Brown et al. using flat, co- planar 
electrodes of different widths and gaps [10| . as shown in 
Fig. [TJ and this general design has been studied exten- 
sively by several groups [H E3, EE EE EE HE Ell- 
in this geometry, the symmetry breaking causes the flow 
to go from left to right by upsetting the delicate bal- 
ance between leftward and rightward slip patterns. The 
fact that pumping comes from a relatively small bias of 
opposing, co-planar surface flows makes this design sus- 
ceptible to flow reversal if nonlinearities in charging dy- 
namics tip the balance the other way, as we demonstrate 
below. Current designs using stepped, three-dimensional 
electrodes are much faster and more robust against flow 
reversal, since the opposing slip patterns work together 
to dri ve p um ping in one direction, as a "fluid conveyor 
belt" [IE EE EE [581 • The sensitivity to flow reversal in 
the planar pump, however, allows us a better opportunity 
to explore subtle nonlinearites in double-layer charging 
dynamics, which is why we have chosen it for our theo- 
retical study. 

In Fig. [5] we demonstrate how the geometrical asym- 
metry results in forward flow when driven at the charging 
time and how steric effects can produce high frequency 
reversal. In Fig. [|3a) and (b) we show the time averaged 
streamlines (a) and time average slip velocity (b) for the 
ACEO pump driven at the charging time, corresponding 
to point A in Fig. 2] We notice a few features of the 
slip profiles. First, the magnitude of velocity over the 
small electrode is greater than the velocity over the large 
electrode. This difference is due to the fact that the total 
charge contained on both electrodes must be the same, 
the local charge density (and thus electroosmotic slip) is 
always higher on the small electrode. Second, we notice 
slip profiles over either electrode is asymmetric with the 
greatest slip found at the edges at the small gap. This 
effect is simply due to high electric field in the small gap 
which can exert a larger force. Finally we see the asym- 
metry in the slip is greater for the large electrode than 
the small electrode, an effect due to the geometry. The 
result for this geometry and frequency is the net flow goes 
in the positive x direction. The relative contributions to 
the flow is easily seen by the cumulative spatial integral 
of the time averaged slip; the dotted line in Fig. [5] (b). 

In Fig. [5] (c) and (d) we show the behavior at the ap- 
proximate peak of the reverse flow in Fig. 0] correspond- 
ing to point B in Fig. [4j The same generic features of 
the slip profiles are the same as in (b). However, the rel- 
ative contributions have changed at this frequency. The 
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FIG. 5: Basic mechanism for flow in ACEO devices. Here we consider a symmetric array of electrodes under a suddenly applied 
voltage. In (a) we show electric field lines at the instance the field is turned on. The electric field points from the positive to 
the negative electrode. The location on the electrodes at y = is shown. In (b) we show the charge density over the electrodes 
at 1/10 the charging time t c , at the charging time, and at 10 times the charging time. In (c) we show streamlines of the electric 
field at the charging time. In (d) we show electric field lines at the charging time. The arrows denote the flow direction. 



reason for the reversal can be explained by the change 
in charging time for the different sized electrodes. When 
steric effects are included, we must remember that the 
capacitance decreases as more charge is stored inside the 
double layer. Since the charge density over the small 
electrode is much greater, the capacitance of the small 
electrode's double layer is decreased at high voltage due 
to the steric effects. The lower capacitance means that 
the charging time of the small electrode is decreased rela- 
tive to the large electrode. This change in charging time 
means that at high frequency and high voltage the small 
electrode has more time to sufficiently charge to generate 
flow. 

This idea is illustrated in Fig. [6](c) and (d). The cumu- 
lative spatial integral of the time-averaged slip shows that 
the negative flow contribution from the small electrode is 
nearly twice the positive contribution. At this high fre- 
quency the large electrode is not very highly charged and 
its forward flow component adds only a small amount to 
the net. That this reversal is caused by the change in 
the charging dynamics can be seen in Fig. [S](e) and (f), 
which shows the same condition as (c) and (d) only with 
the linear Debye-Huckel model; point C in Fig. 0J Here, 
since the large electrode has the same charging time as 
the small electrode, its forward component is significant 
to provide a net positive flow. It is clear from these pic- 
tures that with the inclusion of steric effects, the direction 
and magnitude of flow in ACEO pumps is determined by 



a delicate balance influenced both by the geometry and 
the charging dynamics. 

Up to now, we have taken the volume fraction, v 1 to 
be a constant. To better see the effect of v on the pre- 
dicted flow, in Fig. [7] we show contours of the flow rate 
as a function of voltage and frequency at different values 
of v. For a fixed physical value of a, different values of 
v represent how the flow would change as a function of 
concentration (in dimensionless terms). The heavy lines 
denote the region where the flow is negative, indicating 
a reversal of flow. We see that the frequency response 
of the flow does not change with voltage at low voltage, 
as would be predicted by the linear Debye-Huckel model. 
At low voltage, we see only a single peak in the forward 
flow at the characteristic frequency. As the voltage is in- 
creased the nonlinearity begins to set in which causes the 
peak of the forward flow to move toward lower frequency 
as predicted by dilute solution theory. As the voltage is 
further increased the steric effects begin to become im- 
portant with the reversal eventually setting in. 

The change in frequency of the maximum forward flow 
at higher voltage as v changes is easily understood. As v 
increases, the value of the maximum value of the capaci- 
tance decreases, shifting the peak response to higher fre- 
quencies. For a given v, as voltage continues to increase, 
the continual decrease in the double layer capacitance 
can be seen in the shift of the maximum of the positive 
flow to higher frequency at high voltage. When there are 
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FIG. 6: Streamlines (a,c,e) and slip velocity (b,d,f) for three different cases. The first row is for V = 100, v = 0.01, and ui = 0.7. 
This frequency corresponds to the peak in the forward flow for the steric model. In (a) we show the streamlines. In (b) we 
show the time-averaged slip velocity (solid line) and the cumulative spatial integral of this velocity (the dashed curve). The 
second row, figures (c) and (d) is for V = 100, v — 0.01, and u> — 3. This frequency corresponds to the peak in the reversed 
flow. The third row, figures (e) and (f) is for V = 100 and to = 3 with the linear Debye-Huckel model 9d = —q- 



no steric effects (or Stern layer) the capacitance diverges, 
which can be seen in Figure W[d) . The inclusion of the 
constant Stern layer capacitance causes this decrease to 
limit at the constant Stern layer capacitance [l9j . 

Qualitatively, our model can predict a flow response 
that looks very similar to experimental results, however 
the quantitative ability of our model is still problematic. 
In Fig. [5]Ja) we show the model predicting the experi- 
mental data of Studer et al. [l2|. In Fig. [5] we treat the 
parameter v as a fitting parameter to approximate the 
data of Figure 6 in Studer's experimental work We 
find a value of v = 0.01 provides a reasonable fit to their 
experimental data. However, given a concentration of 0.1 
mM KC1, this value corresponds to a = 4.4 nm, about an 
order of magnitude larger than we would expect based 
on the hydrated radius of the ion. The fact their ex- 
periments [l2| only predicts reverse flow at high voltage 
(their Fig. 7) is consistent with the low frequency regime 
of their work being unaccessible due to electrochemical 
reactions and degradation of the electrode (see their Fig. 
4). 

We considered whether the Carnahan-Starling (CS) 



form of the excess chemical potential would provide a 
more reasonable value of a. In order to include the CS 
equation in the model we numerically solve for the re- 
lationship of (q) to replace Eq. [7] in our model. We 
follow the work of Bieshuvel and Soestbergen [45j to ob- 
tain this relationship numerically, as described in the Ap- 
pendix. The result of using the CS equation in our ACEO 
model is shown in Fig. [8] (b), and it is qualitatively simi- 
lar to our prior results with the Bikerman model. The CS 
steric model has somewhat different quantitative behav- 
ior but still requires an ion size of approximately a = 2.2 
nm to fit the experimental data. While the magnitude of 
the flow reversal changes with the CS equation, the use of 
a more accurate model for the excess chemical potential 
still require a relatively large ion size. 



V. DISCUSSION 

In order to fit the experimental data, the ACEO mod- 
els modified for steric effects imply an ion length scale, 
a, almost an order of magnitude larger than we would 
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FIG. 7: Flow response as a function of voltage for various values of v, the steric parameter. We show values of u = 0.01 (a), 
v = 0.001 (b), v = 0.0001 (c), v — (d), The contours show constant levels of total flow rate. The heavy line shows the U — 
contour and the region above this line denotes the region where flow reversal occurs. 



expect. A factor of 10 in ion size results in a factor of 
1000 in v. This difference means that while the response 
in Fig. \Tia) is closer to what is observed, it is Fig. [7] (c) 
that is predicted based on a reasonable hydrated radius 
for the ions. The current model cannot yet predict the 
experimental data from first principles. It is important 
to note that the large value of v is needed not to predict 
flow reversal at all, but to predict flow reversal at the 
frequencies observed in experiments. 

As discussed in section HH the mean-field models of 
steric effects which we consider are unlikely to accurately 
capture the dynamical formation of a condensed layer on 
the surface. One reason is the difficulty in describing 
the jamming transition, as ions and solvent molecules 
are squashed onto the surface by large, time-dependent 
normal electric fields. Under such extreme conditions, 
we should not expect more than qualitative trends from 
such models. Another possibly crucial feature missing in 
these models is the interaction with the surface, which 
is assumed to be a mathematically flat, homogeneously 
charged, perfectly rigid wall, as sketched in Fig. 2. Since 
the condensed layer of (presumably solvated) ions is at 



the scale of only a few molecular lengths, there must be 
a strong coupling with nanoscale roughness on the sur- 
face. This may have the effect of increasing the effective 
ion size, as we have inferred, since what matters is not 
the ion density in the solution, but the density very close 
to the surface, which is reduced by atomic asperities. 
The viscosity may also diverge near the surface as ions 
become jammed among surface inhomogeneities (29J. It 
would be interesting to do experiments with atomically 
flat electrodes (e.g. carbon graphene) or ones with con- 
trolled nanoroughness to test this hypothesis. 

It is also possible that changes in permittivity, e, within 
the double layer may also effect the frequency response. 
Electrochemists infer for the Stern layer that e may be 
1/10 of e in the bulk [4l[. If it were true that e is 
reduced by 1/10 in the condensed layer, the overall dou- 
ble layer capacitance would be decreased. Such change 
would cause the peaks in the frequency response to shift 
to higher frequency for the same value of v. If we make 
the assumption that e is reduced by a factor of ten in the 
entire double layer, then an ion size of approximately 1 
nm is needed to predict Studer et al.'s [12| data with the 
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FIG. 8: Flow as a function of frequency at different voltages. 
The voltages are V = 0.5, 1, 1.5, 2.9, 3.6, 4.5, and 6 Volts RMS. 
The voltages are selected to match Figure 6 in Studer et al. 
fl2| . The data is presented in dimensional form in order to 
compare the frequency response. The flow velocity is cor- 
rected for the hydraulic resistance of the experimental flow 
loop. In (a), using the Bikerman model, we used a — 4.4 nm. 
In (b), using the Carnahan-Starling model, we used a = 2.2 
nm. The general shape of both frequency response curves 
have many of the experimentally observed features. 



imation breaks down when ion spacings approach the 
Bjerrum length, lg = (ze) 2 / 'AwekT , which is 7 A for 
bulk water and monovalent ions (z = 1). Again, if it 
were true that e is reduced by 1/10 in the condensed 
double layer, then Is = 7 nm. Correlation effects on 
electro-osmotic flow (which to our knowledge have never 
been studied) could be very significant at large voltages, 
even in dilute bulk solutions [291 ]. 

While much work remains, we have highlighted an 
interesting and practically important application where 
dilute solution theory clearly breaks down. We have 
demonstrated that accounting for steric effects in double- 
layer models can have a dramatic impact on the pre- 
dicted flow in ACEO pumps. Because of the large value 
of a needed to fit the experimental data, we cannot claim 
that our model definitively explains high-frequency flow 
reversal in ACEO pumps, but it is the only plausible ex- 
planation to date. The prediction is also quite robust: 
we have found that flow reversal can be predicted by any 
model which accounts steric effects by the generic feature 
that the double layer capacitance decreases with applied 
voltage. 

Physically, it seems that this generic effect of decreas- 
ing capacitance must be true for any model of the double 
layer considering finite sized ions [30( . As the voltage in- 
creases, more ions are packed in the double layer. Once 
the concentration of ions becomes so high that the pack- 
ing density is reached, the ions would have to pile up 
forming a condensed layer whose thickness would increase 
with voltage. As the double layer thickness grows, the ca- 
pacitance must decrease. We believe this general effect is 
true regardless of the details of the model, and we have 
shown that it has major implications for ACEO pump- 
ing and, by extension, all nonlinear ICEO flows. Steric 
effects are clearly among the new effects which must be 
considered in developing more accurate theories of ICEO 
flow at large applied voltages [29j . 



CS equation. If e is reduced only in the condensed layer, 
the true effect on the frequency response would be less 
than predicted with the simple assumption. Therefore, 
it seems likely that accounting for changes in e would re- 
duce the needed ion size to correlate the data, however 
at most this would reduce the ion size to 1 nm, still much 
larger than expected. 

A further complication is that the mean-field approx- 
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Appendix: Modified double-layer models 

In dilute solution theory, the chemical potential of a 
point-like ion i takes the ideal form, 

H i i deal = kT log c t + z t e(j> (10) 

where Cj is the mean concentration and </> is the self- 
consistent mean electrostatic potential, both continuous 
functions in a large dilute system. In concentrated so- 
lutions, the chemical potential is modified by a vari- 
ety of statistical effects, such as electrostatic correlations 
(beyond the mean-field approximation) and interactions 
among discrete, finite-sized ions and solvent molecules. 
In order to isolate such effects, it is customary to decom- 
pose the chemical potential into ideal and excess contri- 
butions, [ii = ^ deal + uf x . Various models for are 
reviewed in Refs. [H M, EH • 
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In the asymptotic limit of thin double layers, the chem- 
ical potential is constant in the normal direction under 
rather general conditions [3!|. In many cases, the alge- 
braic conditions {pt = constant} then suffice to deter- 
mine the concentration profiles {cj} in the diffuse layer 
in terms of the mean potential distribution 4>. For ex- 
ample, in an ideal dilute solution (|10[) . the condition 
Hi =constant yields the Boltzmann distribution Eq. (fT]). 

In the general case, solving the equations {fii = 
constant} and substituting the charge density p(cj>) 
in Poisson's equation leads to a modified Poisson- 
Boltzmann (MPB) equation for the mean potential, 



(11) 



which can be integrated to obtain the differential capac- 
itance of the diffuse layer, 



2s 



la ° Pi<P)d<p 



(12) 



in terms of p((p). This procedure can be carried out 
analytically for the Gouy-Chapman (PB) and Bikerman 
models to obtain Eqs. {2} and (|3|), respectively, or nu- 
merically for more complicated models. 



There are a number of models for the excess chem- 
ical potential due to steric effects of finite (typically 
hydrated) ion sizes. Here, we list a few cases, follow- 
ing Biesheuvel and van Soestbergen [45|]. Bikerman's 
model corresponds to the continuum limit of a lattice 
gas for cubic ions of size a. In terms of the volume fr ac- 
tion <| = y„ . a 3 Ci, the excess chemical potential is sim- 

ply Mm ' 



-fcTln(l - $) 



(13) 



which comes from the entropy solvent molecules (empty 
lattice sites). A more accurate expression for a liquid is 
the Carnahan-Starling (CS) equation [49(, 



$(8-9$ + 3$ 2 ) 

(l-$)3 



(14) 



which has been closely validated by simulations of the 
hard sphere liquid [50(1 • There are also analytical exten- 
sions to mixtures of hard spheres of different sizes, which 
have been applied to volume effects in multicomponent 
electrolytes (45j |. 



